Do Transition Metal Substitutions Dope Carriers in Iron Based Superconductors? 
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We investigate the currently debated issue concerning whether transition metal substitutions 
dope carriers in iron based superconductors. From first-principles calculations of the configuration- 
averaged spectral function of BaFe2 As2 with disordered Co/Zn substitutions of Fe, important doping 
effects are found beyond merely changing the carrier density. While the chemical potential shifts 
suggest doping of a large amount of carriers, a reduction of the coherent carrier density is found due 
to the loss of spectral weight. Therefore, none of the change in the Fermi surface, density of states, or 
charge distribution can be solely used for counting doped coherent carriers, let alone presenting the 
full effects of the disordered substitutions. Our study highlights the necessity of including disorder 
effects in the studies of doped materials in general. 



PACS numbers: 74.70.-b, 71.15.-m, 71.18.+y, 71.23.-k 

Doping is one of the most powerful tools for tuning the 
electronic properties of functional materials. Well known 
examples include n- and p-type semi-conductors, dilute 
magnetic semi-conductors and numerous heavy fermions 
and strongly correlated oxides (e.g. manganites, cobal- 
tates, and cuprates.) In most of these cases, inclusion of 
doped impurities are considered an effective way to intro- 
duce carriers, but recently the amount of carrier doping 
has become the center of a serious debate in the field 
of iron based superconductors (FeSCs). A density func- 
tional theory (DFT) study [1 reported that the addi- 
tional charge density in transition metal (TM) substi- 
tuted FeSCs are fully localized at the TM dopant site, 
which lead [1 to the surprising conclusion that the TM 
substitutions do not dope carriers. While the lack of 
additional density at the Fe sites has recently been con- 
firmed experimentally \3T, '33 , this widely discussed con- 
clusion of the system being undoped, however, seems to 
be in contradiction with several experimental [IHi] and 
theoretical [3 -6 studies that indicate a change of Fermi 
surfaces (FSs) consistent with carrier doping. A serious 
debate was thus initiated in the field: do the TM substi- 
tutions dope carriers into the FeSCs or not? 

In this Letter we investigate this issue via first- 
principles DFT calculation of disordered substitutions. 
Contrary to the previous theoretical studies, which em- 
ployed rather small periodic super cells, we compute the 
configuration-averaged spectral function of disordered 
TM substitutions in BaFe2As2, using the weak Co and 
strong Zn impurities as representative examples. In the 
case of Co substitution, despite the localized charge dis- 
tribution, no apparent localized states are found and a 
large chemical potential shift results from the injection 
of additional carriers to the system. In the case of Zn 
substitution, while both strongly and weakly localized 
states are found, the chemical potential is still shifted 
noticeably, reflecting real carrier doping to the system. 



However, an important loss of coherent spectral weight 
results from scattering against the disordered impurities. 
Consequently, the actual number of coherent carriers in 
the system reduces noticeably from the Luttinger count 
of the FSs. Our findings illustrate that none of the FSs, 
density of states, or the charge density can be used for a 
physically meaningful counting of coherent carriers in a 
disordered system in general. More importantly, physical 
effects of the disordered impurities are much richer than 
simply adding carriers. Specifically, the emergence of in- 
coherent carriers and loss of coherent carriers holds in- 
teresting potentials to strengthen superconductivity over 
magnetism, and should be taken into account in future 
studies. 

The configuration-averaged Wannier based spectral 
function {An{k^uj)) of Wannier orbital n, frequency cj, 
and crystal momentum k (in the 2-Fe Brillouin zone [7]) 
is obtained by averaging over 10 large random-shaped 
supercells containing 400 atoms on average, including 
12.5% random TM substitutions (e.g. Fig. [ija)). (For 
the spectral functions at fixed k-points in Fig. [2|c)-(f) 
and [3]^b)(c) 100 supercells of 800 atoms on average are 
used). The spectral functions are calculated directly 
from the supercell eigenvalues and eigenvectors, by ap- 
plying the recently developed unfolding method [8 . The 
first principles simulations of the large supercells, nec- 
essary for a proper treatment of disorder, are made af- 
fordable by the recently developed Wannier function [10] 
based effective Hamiltonian method for disordered sys- 
tems [TT]. The influence of the TM substitutions on the 
Hamiltonian is extracted from three DFT p!2HT4j calcu- 
lations: the undoped BaFe2As2 and the impurity super- 
cells Ba2Fe3CoAs4 and BagFeisZnAsie- The low energy 
Hilbert space is taken within [-10,3] eV consisting of Wan- 
nier orbitals of Fe-d, Co-d/Zn-d, and As-p characters. 
The accuracy of the calculated effective Hamiltonian is 
benchmarked [14 against full DFT calculations. 
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FIG. 1: (color online) Two seemingly contradicting results: 
charge localization and carrier doping, (a) An example of 
large sized supercells used for the configurational average, 
(b) The average occupation of the atomic Wannier orbitals 
showing the doped charge to be localized at the Co site. The 
FSs of (c) BaFe2As2 and (d) Ba(Feo.875Coo.i25)2As2 showing 
the hole pockets at F to shrink and the electron pockets at 
X to grow, consistent with adding 0.125 electrons/Fe to the 
undoped system. 



Let us first re-examine the previously reported charge 
localization [1 and chemical potential shift [5, 6 within 
our disordered Co substitution case. Fig. [TJb) shows the 
average occupation of each atom calculated from fre- 
quency/momentum integration and orbital summation 
of {An{k^uj)). Comparing the occupations in the doped 
system with that in the undoped system shows that, in 
agreement with Ref. [H [3 |6], the additional charge is 
indeed distributed at the Co impurity, seemingly imply- 
ing that Co does not dope carriers. On the other hand, 
compared with the FSs of the undoped (Fig. [ijc)) sys- 
tem, the doped system (Fig.jljd)) contains smaller/larger 
hole/electron pockets around the F/X point, in agree- 
ment with the previous experimental and theoretical find- 
ings [2 - 6 . Furthermore the FSs of the doped system are 
found to be consistent with the undoped system upon 
adding the nominal doping of +0.125 electrons per Fe 
(the dotted line in Fig. [ijd).) The same can also be 
observed from the band structure in Fig. |2|b). So, our 
results of disordered Co substitutions confirm both the 
localized distribution of additional charge [ij 13 |^ and 
the chemical potential shift [5l [6] , which still leaves us 
with the question: Do TM substitutions dope carriers or 
not? 

To reconcile these seemingly contradicting results, a 
notion of impurity potential Ej (or self-energy ^ ) in- 
duced "effective doping" was previously proposed [5l|6l. 
However, our results give a different physical explana- 
tion and show that such notion is not only unnecessary 
but also negligible in effect for the case of Co substitu- 
tion. Let's focus on the (/c,a;)-dependence of {An{k^uj)) 
in Fig. [2] Notice first that while the band dispersion 
resembles well the undoped system, important effects 
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FIG. 2: (color online) Configuration averaged spectral func- 
tion of disordered Ba(Feo.875Coo.i25)2As2 in the full low en- 
ergy Hilbert space (a) where the partial DOS of Co is en- 
hanced by a factor of 7.(b) Comparison of spectral func- 
tion and DOS of disordered system (black) with the undoped 
one (grey) near the chemical potential. Bloch based spec- 
tral function of band j — 12 at (c) ko — 0.15FM and (d) 
ko = 0.25FM and of bands j = 8, 11 at (e) /co = F and (f) 
/co = 0.5FM, where the dashed/solid lines indicate the cen- 
ters of the dispersive/non-dispersive components. The black 
curves in (c)(d)(e)(f) correspond to the total spectral func- 
tion. 



of disorder scattering can be observed. For example, 
the Ye-d bands (even at the chemical potential) develop 
line widths in both momentum and frequency, reflecting 
their finite mean free path and lifetime from scattering 
against Co impurities. Also notice that besides these 
features, no apparent localized state (flat band) is found 
in Fig [2] to host one electron per Co. Therefore, the 
localized distribution of additional charge must originate 
from slight enhancement of the wave function near the Co 
sites in a large number of extended occupied states, in- 
stead of a small number of strongly localized or resonant 
states! This deduction is further supported by the DOS 
in Fig.[2|a), in which the broad Co spectral weight (x7 for 
better comparison) is shown to be transferred to slightly 
lower energy, in relation to the Fe spectral weight. That 
is, contrary to the previous interpretations [U [3 |6] our 
result shows that practically no carriers are "trapped" 
around Co, but instead a large number of itinerant car- 
riers conspire to enhance the charge distribution around 
Co in order to screen the additional proton. Therefore, 
the notion of "effective doping" is unnecessary in this 
case, and in fact negligible in effect. Since the change 
of the wave functions and the band dispersion are rather 
small in this simplest case, from any physically mean- 



3 



ingful consideration concerning doping, Co substitution 
really just adds about one more itinerant electron per Co 
to the system. 

However, these itinerant carriers are not all the same in 
nature, as some of these carriers lose their coherence due 
to scattering against disordered impurities. The lack of a 
well-defined dispersion indicates diffusive propagation of 
these incoherent carriers, in contrast to the particle-like 
propagation of the coherent carriers. A clear example 
of non-dispersive spectral weight can be found in band 
12 around ko = F. As illustrated in Fig. |2|c)(d), the 
band consists of two components: a peak (dashed line) 
that disperses from 0.6 to 0.7eV along FM and a non- 
dispersive satellite (solid line) that remains fixed a little 
above 0.8eV. Other incoherent features are more difficult 
to recognize due to the presence of multiple bands in the 
same frequency range. To resolve the different bands, in 
Fig. [2|c)-(f) and [3lb)(c), we analyze the spectral func- 
tion Aj{k^uj) = -^lm{kj\G{uj)\kj) on the basis of the 
undoped Bloch states \kj) of crystal momentum k and 
band index j [14]. In Fig. [2|e)(f) an incoherent compo- 
nent in the hole pocket is illustrated. While along FM 
the main peak in band 11 disperses from 0.1 to -0.3eV, 
the satellite around -0.6eV remains fixed. Note also that 
while nearby in frequency, the non-dispersive satellite of 
band 11 is not part of the dispersing band 8. Due to the 
proximity of the incoherent features to the main peaks, 
an unambiguous quantification of the obviously reduced 
coherent spectral weight is challenging in this case. 

From the perspective of disordered impurities, Co sub- 
stitution of Fe introduces only a weak impurity potential 
~-0.5eV to the d-orbitals, much smaller than the band- 
width of the d-bands, and thus can be considered as a 
small perturbation to the system. This is ultimately the 
reason why the Fermi surfaces resemble very much the 
undoped system added with nominal doping of +0.125 
electrons per Fe (e.g. Fig. 1(d)), even though Luttinger 
theorem should not apply to disordered systems. There- 
fore, it is instructive to explore the Zn substitution of Fe, 
which represents the other limit with very strong impu- 
rity potential (~-7eV) much larger than the bandwidth 
of the d-bands. 

As expected, the strong impurity potential of Zn sub- 
stitution is able to induce real localization, in great con- 
trast to the Co substitution. Fig. |3|a) shows our result- 
ing band structure of Zn substituted system with fiat 
core levels of Zn-d character, which are good examples 
of strongly localized states created by the large impurity 
potential. More interestingly, weakly localized As-states 
can also be observed clearly within [-6,-5] eV, right be- 
low the As-p bands [31]. The small region of /c-space in 
which these weakly localized states reside reflects their 
large size compared to the 2-Fe unit cell (similar to the 
oxygen vacancy states in ZnO [23]). 

Due to the very strong impurity potential, the in- 
fluence of the impurity can no longer be considered as 
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FIG. 3: (color online) Configuration averaged spectral func- 
tion of disordered Ba(Feo.875Zno.i25)2As2 in the low energy 
Hilbert space (a). Bloch based spectral function (b) of band 
j = 11 at /co = r and (c) of band j = 12 at ko = 2/3MX. 
(d) FSs of doped system compared with the undoped system 
with ^ 0.2 additional electrons per Fe. (e) Band structure 
and DOS of doped (black) and undoped (grey) system near 
the chemical potential, (f) The estimated number of doped 
coherent carriers. The black curves in (b)(c) correspond to 
the total spectral function. 



a small perturbation, and consequently the counting 
of the carriers would become non-trivial. Indeed, the 
broadening of the Fe-d spectrum in Fig. [sj^b) is signifi- 
cantly stronger than the previous case, indicating a much 
shorter lifetime and mean free path. In fact, even As-p 
bands scatter strongly in Fig. [sja). Now, from a com- 
parison of the valence of Fe^+ and Zn^+, one could have 
expected that Zn substitutions do not dope carriers to the 
system. This is obviously not the case, since the chemical 
potential has been shifted so much that the hole carriers 
are almost entirely removed from the system. The new 
FSs in Fig. |3|d) appear to have been doped with 0.206 
electrons per Fe, a number that cannot be easily ratio- 
nalized from simple counting of valence, in great contrast 
to the case of Co substitution. This illustrates clearly the 
general inapplicability of Luttinger theorem to disordered 
systems. 

Of course, the physically meaningful count should dis- 
tinguish the coherent carriers from the incoherent fea- 
tures. For this purpose we analyze the components of 
the Bloch based spectral function. Two examples are il- 
lustrated in Fig. |3jb)(c), representing the bands in the 
hole and electron pockets respectively. An integration of 
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the Bloch based spectral function within [-l,l]eV gives 
an estimated coherence factor of ~0.85. This results in 
an effective doping of ~+0.049 / '^-0.126 per Fe in the 
electron/hole pockets (see Fig. [sjf)). A simple Luttinger 
count would lead to a large overestimation of 73% of the 
number of doped electron carriers. 

Interestingly, the Co/Zn substitutions are more effi- 
cient in removing the coherent hole carriers than adding 
the electron ones, an effect relevant to the doping depen- 
dence of Hall coefficient p4H26] . Of the ^0.15 spectral 
weight loss with Zn substitutions, ~ 0.09 comes from the 
trivial strongly localized Zn d-orbitals, while the remain- 
ing ~ 0.06 corresponds to the decoherence of Fe bands. 
(Similar spectral weight losses can be expected in the case 
of Fe vacancies.) Finally, Fig. |3|e) shows an enhanced 
DOS upon substitution, even though the coherent frac- 
tion of DOS actually decreases. This is quite alarming, 
considering that DOS is routinely used in the standard 
studies of doped materials, but actually provides no dis- 
tinguishability of coherent features from the incoherent 
ones. 

The above disorder effects have important physical im- 
plications that immediately impact on experimental and 
theoretical studies of FeSCs. At the simplest level, one 
expects a doping dependence of the broadening of the 
bands, as found in the raw ARPES data \T, ^3 , and 
contribution to the incoherent features in other spec- 
troscopies [151417] . But, deeper issues are present. As 
demonstrated above, other than the appearance of inco- 
herent features, the chemical potential and the associated 
phase space of coherent carriers are no longer related to 
the coherent carrier density. Indeed, we have found the 
Fermi velocities of the electron/hole carriers to resemble 
those of a highly doped system, but instead with a much 
smaller coherent carrier density. This cannot be possibly 
described with merely adding carriers to the undoped 
system, and is in contradiction with the textbook theo- 
ries, based on which most experiments are analyzed. For 
instance, in optical conductivity and dielectric function, 
one might now find a spectral line shape similar to that 
of an undoped system with additional electrons (since its 
band dispersion resembles the real disordered system), 
but with a overall reduction in the intensity (due to the 
weight loss of coherent carriers) as dictated by various 
sum rules. The same consideration should apply to trans- 
port and other spectroscopies as well, since they are all 
sensitive to both the phase space controlled by the band 
structure and the spectral weight related to the coherent 
carrier density. 

More importantly, besides the impurities' standard 
role of suppressing Tc via scattering [27H29] , the emer- 
gence of incoherent carriers and the lost of coherent spec- 
tral weight can also affect the competition between mag- 
netic and superconducting instabilities in the underdoped 
regime. Both effects serve to suppress long-range antifer- 
romagnetic order, allowing a stronger superconductivity 



to develop. This is because the the nesting condition 
is reduced by the loss of coherent spectral weight, while 
the ill-defined momentum of the incoherent carriers give 
negligible contribution. This explains, at least in part, 
why both electron doping and hole doping both suppress 
the antiferromagnetic order, despite the degree of nesting 
condition. On the other hand, the zero momentum su- 
perconductivity still benefit from all the incoherent car- 
riers. In fact, the slow propagation of the incoherent car- 
riers makes them more susceptible to the strong coupling 
regime of superconductivity [30 that hosts stronger local 
pairing but weaker phase coherence. (This might offer a 
resolution to the current issue of power law behavior of 
the penetration depth [18 ). Coupled to the lighter co- 
herent carriers, this can actually help enhance the super- 
conductivity by providing a stronger local pairing. Such 
opposite effects on the magnetism and superconductivity 
have in fact been observed recently in Zn substitution of 
F-underdoped LaOFeAs [26j. Finally, the loss of coher- 
ent spectral weight would allow practically the current 
itinerant pairing theories p!9H22 ^ to incorporate stronger 
interaction parameters without triggering a long-range 
magnetic order, and thereby to further enhance the su- 
perconducting instability. 

In conclusion, by computing the configurational- 
averaged spectral function of disordered Co/Zn substi- 
tuted BaFe2As2, we identify important physical effects 
of TM substitution that cannot be captured by merely 
adding carriers to the undoped systems. While both 
strongly and weakly localized states are found via Zn 
substitution, no apparent localization is found via Co 
substitution, despite the localized distribution of addi- 
tional charge. On the one hand, the FSs and Fermi ve- 
locities are modified as a consequence of chemical po- 
tential shift, reflecting injection of doped carriers. On 
the other hand, some incoherent carriers emerge, in as- 
sociation with a loss of the coherent spectral weight that 
reduces considerably the actual number of coherent car- 
riers. Our results illustrate clearly that generally in dis- 
ordered systems, none of the FSs, DOS, or the charge 
distribution can be solely used for the counting of co- 
herent carriers. Physically, our reported features would 
produce unusual behaviors in various spectroscopies and 
other physical properties beyond the standard textbook 
descriptions, and thus needs to be taken into account 
in future experimental and theoretical studies. In addi- 
tion, they offer new possibilities to tip the balance from 
magnetism toward superconductivity, an important as- 
pect that deserves to be considered in future studies of 
superconductivity in general. 

Work funded by the U S Department of Energy, Of- 
fice of Basic Energy Sciences DE-AC02-98CH10886 and 
DOE-CMCSN. 
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Supplementary information 



DETAILS OF THE DENSITY FUNCTIONAL THEORY CALCULATIONS 
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FIG. 4: (Left) undoped normal cell BaFe2Co2, (middle) single Co impurity supercell Ba2Fe3CoAs4 and (right) single Zn 
impurity supercell BasFeisZnAsie- 



The lattice structure of the undoped normal cell, BaFe2As2, is depicted in the left of figure |4] Its space group 
I4/mmm and its lattice parameters a=7.49 bohr and c=24.60 bohr and z=0.3545, were taken from [2 . The primitive 
normal cell lattice vectors expressed in Cartesian coordinates are given by: ai = —^ax-\-^ay-\-^cz ; a2 = ^ax — 
^ay + ^cz ; as = ^ax + ^ay — ^cz. The supercell used to capture the influence of the Co impurity Ba2Fe3CoAs4 is 
depicted in the middle of figure [4] and its primitive lattice vectors expressed in normal cell lattice vectors are given by: 
Ai = a2 -\- as '■, A2 = ai -\- as '■, A3 = ai + a2. To capture the influence of the Zn impurity a four times larger supercell 
is used to properly treat its significant non-local influence on the nearest As and Fe orbitals. The BagFeisZnAsie 
supercell is depicted on the right of figure |4] and its primitive lattice vectors expressed in normal cell lattice vectors 
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are given by: Ai = 2a2 + 2as ; A2 = 2ai + 2as ; A3 = ai + a2. We applied the WIEN2K[T] implementation of the full 
potential linearized augmented plane wave method in the local density approximation. The k-point mesh was taken 
to be 10x10x10 for the undoped normal cell and 10x10x3 and 5x5x3 for the Co and Zn supercells respectively. 
The basis set sizes were determined by RKmax=7. 



BENCHMARKS OF THE EFFECTIVE HAMILTONIAN AGAINST DFT 



To explore the accuracy and efficiency of the effective Hamiltonian for the case of Ba(Fei_a;Ma;)2 As2 with M=Co/Zn, 
we present comparisons of spectral functions An{k^uj) calculated from the full DFT and the effective Hamiltonian 
(see figures |5][7| . The size of the deviations between the full DFT and the effective Hamiltonian should be compared 
with the size of the impurity induced changes. For this purpose the spectral function of the undoped BaFe2As2 
is also plotted as a reference for each benchmark. The test cases in figures |5] and [6] are designed as extreme test 
cases for the "linearity" approximation of the impurity influence (see formula (1) of Ref. [3 ). The case of BaCo2As2 
corresponds to the maximum extrapolation, from x = 0.25 in the single impurity cell (see middle of Fig. [4| to x = 1. 
For the Zn doping, the BaFeZnAs2 system is more meaningful case then BaZn2As2, in order to have remaining Fe 
bands. In this benchmark the dispersion of the Zn bands around -8 eV are less satisfactory captured by the effective 
Hamiltonian. This can be understood as a good example in which additional two-impurity effects are needed in Eq.(l) 
of Ref. [3 . Nonetheless, the more relevant Fe bands near the chemical potential are very reasonably described, despite 
the extreme extrapolation in doping from x = 0.0625 in the single impurity cell (see right of Fig. [4| to x = 0.5. The 
test case in figure [t] of a a/5 x supercell is designed to test the partitioning of the impurity influence from its super 
images (see section IV of the supplementary of Ref. [3 ). The basis set of Linear Augmented Plane Waves (LAPW's) 
used in the full DFT is ^ 30 times larger then the basis set of Wannier functions used in the effective Hamiltonian 
method. Since the number of floating point operations of diagonalization depends cubically on the size of the matrix 
this implies an efficiency increase by a factor of 30^ ^3-10^. Furthermore the full DFT calculation involves multiple 
self consistent cycles (typically 15) whereas the effective Hamiltonian method (as currently implemented) requires 
only a single diagonalization, which increases the efficiency by another order of magnitude to ~ 6 • 10^. 
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CONVERGENCE WITH RESPECT TO SIZE AND NUMBER OF CONFIGURATIONS 



In figures |8]/|9] we demonstrate the convergence of the configurational averaged spectral function {A{k = ko^uj)) of 
Ba(Feo.875Zno. 125)2 As2 at ko = T/ko = |MX, corresponding to the hole pocket/weakly localized As states. 




FIG. 8: Demonstration of convergence of the configurational averaged spectral function of Ba(Feo.875Zno.i25)2As2 for the 
energy distribution curve {A{k = ko^uj)) at ko = P. 
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FIG. 9: Demonstration of convergence of the configurational averaged spectral function of Ba(Feo.875Zno.i25)2As2 for the 
energy distribution curve {A{k ko,u)) at A;o |MX. 



BLOCK BASED SPECTRAL FUNCTION 



The Bloch based spectral function is obtained from: 

A^{k,uj) = --lm{k3\G{uj)\kj) = Y^5{uj- e^j)|(fci|/f J)|2 ; {kj\KJ) = Y^Mkn) {kn\K J) (1) 



with k/K is the normal/super ceh crystal momentum, cu the frequency, j/J the normal/super cell band index, n the 
normal cell Wannier index. Here ckj and \KJ) are the supercell eigenvalue and eigenstate, {kn\kj) is the Wannier 
based Bloch state of the undoped system and {kn\KJ) is the unfolding matrix element (see formula (4) of Ref. [4|). 



In figure 10 we display the Bloch based spectral function for Ba(Feo.875Zno.i25)2As2 at ko = |MX. 
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FIG. 10: Bloch based spectral function {Aj{k 
I MX 



ko 



ko,u)) of band index j for Ba(Feo.875Zno.i25)2As2 at the fixed k-point 
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